2D transport and screening in topological insulator surface states 
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We study disorder effects on the surface states of the topological insulator Bi2Se3 close to the 
topologically protected crossing point. Close to charge neutrality, local fluctuations in carrier density 
arising from the random charged disorder in the environment result in electron and hole puddles that 
dominate the electronic properties of these materials. By calculating the polarizability of the surface 
state using the random phase approximation, and determining the characteristics of puddles using 
the self-consistent approximation, we find that band asymmetry plays a crucial role in determin- 
ing experimentally measured quantities including the conductivity and the puddle autocorrelation 
length. 
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Topological insulators (TIs) are a new class of materi- 
als [lfy that are typically distinguished by their robust 
metallic surface state encapsulating a non-conducting 
bulk. From this perspective, a major shortcoming of the 
early experiments was that carrier doping levels were suf- 
ficiently high, resulting in the bulk bands being metallic, 
as opposed to insulating. This made it almost impossible 
to definitively separate the properties of the surface state 
from that of the bulk, and the system did not behave as 
a TI because the bulk is conducting rather than insulat- 
ing. This leads to the problem that while spectroscopic 
measurements such as ARPES demonstrate the clear ex- 
istence of the expected surface topological bands, trans- 
port measurements have been difficult to interpret due 
to the coexistence of both bulk and surface conduction. 

It is, therefore, encouraging that several recent exper- 
imental studies [5|48| report the direct observation of 2D 
surface states in transport measurement. These reports 
claim to observe the electronic properties of the surface 
state with energy close to the topologically protected 
band crossing point (also called the Dirac point). Our 
previous work on the electronic properties of graphene [9( 
would lead us to expect that close to the Dirac point 
the energy landscape and the spatial electronic structure 
would become highly inhomogeneous, breaking the sur- 
face into puddles of electrons and holes. The charge in- 
homogeneity would also result in a low density plateau 
in the conductivity with a non-vanishing minimum con- 
ductivity at the Dirac point [§- 12 1. Given the highly dis- 
ordered nature of the experimental systems, we expect 
the electron-hole puddles to completely dominate the 2D 
transport and electronic structure near the Dirac point. 

This physically intuitive picture of inhomogeneous 
electron and hole puddles has been highly successful in 
understanding both transport experiments 13|, ll 1| and 
scanning probe experiments [lEj in graphene. The goal of 
this work is to understand the role of a locally fluctuating 
carrier density on the electronic screening and conductiv- 
ity of the 2D surface states in these new 3D TI materials. 
Our main finding is that band asymmetry, which can al- 



most always be ignored when studying graphene, plays 
a crucial role in TI experiments. One striking example 
is that in some cases, the TI surface conductivity does 
not even have a minimum conductivity in the vicinity of 
the band-crossing point, making the physics very differ- 
ent from that of graphene. Comparing our theory with 
the existing Bi 2 Se3 experiments, we also conclude 

that the current systems have very large background dis- 
order making it difficult to access the physics of the Dirac 
point. 

The starting point for our calculation is to make some 
reasonable approximation for the band structure of the 
topological surface state. One approach 1|| 17 1 would 
be to perform ab initio calculations and fit that data to 
the most general model Hamiltonian allowed by symme- 
try However, we find significant discrepancy between the 
electronic structure calculations and the photoemmision 
(ARPES) experiments. As a result, a reasonable com- 
parison with experiment would require us fine-tuning a 
model Hamiltonian with 12 parameters. While this pro- 
cedure could be done, it would be unnecessarily cumber- 
some and obscure any physical insight. Instead, we follow 
the minimal model proposed earlier 18j in the literature 
mimicking the full ab initio band structure 
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where (o~ x ,o~ v ) is a 2D vector of Pauli matrices, k — 
(k x ,k y ) is the 2D wave vector, vf is the Fermi velocity 
of the Dirac bands, and effective mass m* characterizes 
the degree of asymmetry between the electron and hole 
bands. Estimates for the values of these two parameters 
for Bi2Se3 vary widely in the experimental literature. For 
example, values for vf vary from 2 x 10 5 m/s [l9[ to 
6.4 x 10 5 m/s 0], and measured values for m* vary from 
0.11 m e [20] (m e is the electron mass) to 0.32 m e plj ]. 
This situation should be contrasted with graphene, where 
vp is the single band parameter, and most experimental 
reports agree on its value to within 5 percent Q . We use 
the two parameter (m* and vp) model of Eq. ([T]) in the 
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FIG. 1. (Color online) Theoretical calculations for dependence of conductivity on carrier density using parameters r s = 0.1 
and no — 2 x 10 13 cm -2 . The left panel has impurity concentration ni mp = 10 13 cm -2 , and the right panel is dirtier with 
iiimp = 5 x 10 13 cm -2 . Dashed blue lines show the Boltzmann transport result described in Eq. (Q, and dash-dotted black 
lines are for Dirac bands i.e. m* — s> oo in Eq. (TTJ). The red solid line is the result of the effective medium theory calculation 
obtained by solving Eq. ((6| using d — 0.1 nm. The left panel has a conductivity minimum in the vicinity of the Dirac point, 
while for the right panel the conductivity monotonically increases as the carriers are tuned from holes (negative n) to electrons 
(positive n). 



current work. 

Since we are concerned with the screening properties of 
electrons, it is useful to define an interaction parameter 
r s = e 2 /(kHvf), where we reiterate that throughout this 
work, v-p is the parameter in Eq. ([T]) characterizing the 
Dirac-like bands; only at low carrier density does it co- 
incide with the Fermi velocity. Here k is approximately 
half the dielectric constant of the bulk E^Sea insulator 
(whose reported value varies from around 30 [19( to about 
55 [12] )■ We certainly expect that as more experiments 
on TI materials become available, these parameters will 
become better known both for Bi2Se3 and other related 
materials. To make our theory more compact, we formu- 
late everything in terms of r s and a characteristic density 
tiq = (m* vp) 2 1 (AttH 2 ) , where reasonable values of r s are 
in the range 0.05 to 0.5, and reasonable values of no are 
from 10 11 cm -2 to 3 x 10 13 cm -2 . We note that tiq is an 
important parameter characterizing the deviation of the 
system (for n > uq) from purely Dirac-like behavior — 
in graphene, no is very large. 

The Thomas- Fermi screening theory for electrons spec- 
ifies that all external potentials are screened by a surface 
2D dielectric function e(q) — 1 + <7tf/<7- For the Hamil- 
tonian in Eq. (JT|), we find 
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where sgn is the signum function, and we use the conven- 
tion that electrons have sgn(n) = +1 while holes have 
sgn(n) = — 1. For both electrons and holes, we have 
fcp = \/|47rn[. Note that <7tf diverges for n — > —no, 
implying perfect screening associated with the diverging 



density of states in Eq. ([T]) arising from the quadratic 
dispersion and the band asymmetry. It turns out that 
the theory can be completely characterized by the two 
parameters r s and no, rather than the three microscopic 
parameters m*, vf, and n. 

Our numerical analysis using the full random phase ap- 
proximation (RPA) shows that the Thomas-Fermi anal- 
ysis is accurate provided we restrict the carrier den- 
sity for holes to |n| < |no|. No such restriction is re- 
quired for electrons. Within the Boltzmann transport ap- 
proximation, the conductivity a = (fcp^)(e 2 /2/i), where 
I = vft/t) is the mean- free path. The scattering time is 
calculated within Born approximation as 
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where 



"amp 



is the surface density of random charged im- 



purities, e k is the carrier energy, kk < is the scattering 
angle between wave vectors k and k', and v(q)/e(q) is 
the Fourier transform of the screened impurity potential. 
For the purpose of this calculation we assume that the 
dominant scatterers are long-ranged Coulomb impurities 
(with an average 2D density of n lmp placed an average 
distance of d away from the TI surface) although our 
formalism can easily be generalized to other types of im- 
purities. For these charged impurities, the conductivity 
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can be calculated analytically, giving 
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The results for the Boltzmann transport theory are 
shown in Fig. [T] One immediately observes that the 
asymmetry in the conductivity is quite pronounced com- 
pared to a linear Dirac dispersion (m* = oo; also shown). 
The electron branch has a much larger conductivity with 
tTs(n) being super- linear, while the hole branch has a 
much lower sub-linear as^n). This pronounced asym- 
metry between electron and hole transport, following di- 
rectly from the band asymmetry of Eq. ([!)). is a charac- 
teristic feature of 2D TI transport. 

At low carrier density, the disorder induced fluctua- 
tions in carrier density become larger than the average 
carrier density. In particular, when the average carrier 
density vanishes with the chemical potential at the Dirac 
point, one might expect that the electronic properties 
of the system are determined by the typical carrier den- 
sity inside the electron and hole puddles. For example, 
one could define a carrier density distribution function 
P[n], and the condition of zero average carrier density 
implies the vanishing of the first moment of P[n\. The 
second moment of the carrier density distribution n rms 
would then determine how the carriers in this inhomo- 
geneous system screen any external potential, including 
the impurity potential that induced the fluctuations to 
begin with. This implies that the density fluctuations 
need to be calculated self-consistently[ll[. We calculate 
the properties of this inhomogeneous system by assum- 
ing a global screening function that depends on the im- 
purity profile only through an effective carrier density 
n e ff (where knowing n e ff, one can then calculate all mo- 
ments of P[n] including n rms , e.g. for Dirac fermions, 
« V3n cS ). 

This effective carrier density n c g is nothing other than 
a measure of the typical carrier density inside the elec- 
tron and hole puddles. After obtaining n e s , we can then 
compute other properties of the Dirac point including its 
conductivity (that can be measured in a transport exper- 
iment) or its density-density correlation function (mea- 
sured in STM). Calculating n c s is therefore a central 
result of this work. We do this by requiring that the 
density induced by the second moment of the screened 
disorder potential is precisely the same as the density 
entering the global screening function. Applying this pro- 
cedure to Eq. (fTJ), and defining the dimensionless variable 
y = |n c ff|/no, we derive a system of equations that can 
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FIG. 2. (Color online) Theoretical calculations for disorder 
induced carrier density n c ff as a function of impurity concen- 
tration nimp using parameters r a = 0.1, no = 10 13 cm -2 , and 
d = 0.1 nm. The curves were calculated using the Thomas- 
Fermi screening theory, while the circles (holes) and squares 
(electrons) show the RPA. For electrons, the two approxima- 
tions agree to within our numerical accuracy. For holes, they 
agree only for n c ff < no (shown as a dotted line). One im- 
plication of this result is that charge neutrality point will be 
distinct from the Dirac point. 



be easily solved numerically 
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(s = — 1) bands. 

In Fig. [5] we show the effective carrier density obtained 
using this self-consistent method. Notice that the in- 
duced carrier density is different for electrons and for 
holes except in the limiting case ni mp <C no [where Eq. ([1]) 
gives a symmetric linear Hamiltonian] . This asymme- 
try implies that the charge neutrality point (where num- 
ber of electrons and holes are equal) does not necessar- 
ily coincide with the Dirac point (the crossing point be- 
tween electron and hole bands). The numerical value of 
n e ff/ni mp is rather small (less than 0.03). This not only 
guarantees convergence of the theory, but also has im- 
portant implications for experiments. Since no and ni mp 
are comparable in current experiments, if n e s /ni mp is not 
small, this would imply that the energy scale associated 
with the disorder-induced inhomogeneity would be com- 
parable to the bulk band-gap, and there would be no hope 
of observing any physics associated with the topological 
surface state. We note that in graphene n c s ~ ni mp 
in contrast to 2D TI transport where we are finding 
n G ff < n imp . 
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FIG. 3. (Color online) Minimum conductivity as a function of system parameters (a) impurity concentration ni mp , (b) dielectric 
constant k, (c) low energy Fermi velocity vf, and (d) impurity distance d. The values for the parameters not being varied are: 
ni mp = 5 x 10 13 cm" 2 , vf = 4.5 x 10 5 m/s, « = 50, d — 0.1 nm, and m* = 0.2. A zero value for tr m in means that there was no 
local minimum in the conductivity in the vicinity of the Dirac point. Solid curves are a guide to the eye. 



We now proceed to calculate the electronic transport 
in this inhomogeneous puddle-dominated carrier density 
landscape. The procedure is to first assume that the lo- 
cal conductivity can be calculated from the local carrier 
density using Eq. ((4]). Then, an effective medium theory 
(EMT) can be used to calculate the global conductivity 
from this distribution of spatially fluctuating local con- 
ductivities [HI 

f j or l ®b(ji) cemt „ , PA 
dnP[n] — — =0. (6) 

J er B (n) + (Temt 

The EMT assumes that the dominant contribution to the 
resistivity arises from scattering inside the electron and 
hole puddles, and not across the puddles. This is a rea- 
sonable assumption because the cross puddle backscat- 
tering is suppressed due to spin conservation. To sim- 
plify the theory we assume that P[n] is Gaussian, with 

variance \/3n* , where n* = »J n^ s n^ s is the geometric 
mean of the electron and hole density fluctuations (see 
Fig. [5]). Our results are shown in Fig. [T] An important 
result of performing the EMT average using an asym- 
metric conductivity is that the carrier density at which 
the conductivity is minimum is distinct from both the 
charge neutrality point and Dirac point, further com- 
plicating the analysis and interpretation of experiments. 
This is illustrated dramatically in the right-hand panel, 
where there is no conductivity minimum in the vicinity 
of the Dirac point. Rather, the conductivity increases 
monotonically from a = 0, for n = — rt to tTs(|n|) for 
n > hq. Indeed, in the experiments of Ref. H, the most 
disordered sample does not show a (clear) minimum con- 
ductivity, while other devices do have a local cr m i n close 
to the Dirac point, similar to that shown in the left-hand 
panel of Fig. [TJ The possible non-existence of a conduc- 
tivity minimum in 2D TI transport is a new qualitative 
prediction of our theory. 

We have explored how tT m i n depends on a variety of 
physical parameters. In most cases, we find that the 
minimum conductivity increases with increasing k, v-p, 



4 

C(r) 3 



2 



\ — n = 







\ e 






\\ — n = 


10 13 


cm 2 


n\ e 






\\\ — n = 


10 15 


-2 

cm . 


\\ x e 









5 10 15 20 

r [nm] 



FIG. 4. (Color online) Theoretical calculations for potential 
correlation fluctuation C(r) = {(V(r) - V)(V(0) - V)) for 
the electron band using the same parameters as in Fig. [lb. 
The y-axis is normalized by (r s Hvf ^/7rni mp ) 2 . From top to 
bottom, the curves represent average electron doping of n — 0, 
10 13 cm" 2 , and 10 15 cm -2 . 

to* and sample purity (n^ p ). However, there is also a 
large range of parameter space where there is no con- 
ductivity minimum associated with the Dirac point (we 
denote this = 0). We illustrate our results 

in Fig. [3l where we show the dependence of (T m ; n on k, 
vf, and d. As seen in the figure, in each of these cases, 
the crossover between the presence or absence of a well 
defined cr m i n occurs quite sharply as a function of the 
parameter being varied. 

The analytic results provided in Eq. ([5]) were done us- 
ing the Thomas-Fermi (TF) approximation. We have 
also done calculations using the full static RPA. We find 
in general, as shown in Fig. 2, very good agreement 
between the TF and RPA results with the RPA results 
being obtained completely numerically. 

Finally, we turn our attention to the recent exper- 
iments which have directly observed the electron-hole 
puddles on the Bi 2 Se3 surface [3]. In addition to be- 
ing a local probe, by tuning the bias voltage between 
the tip and the sample, these experiments can probe fea- 
tures in the density of states away from the Fermi energy. 
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In current TI materials, this is especially important be- 
cause unintentional doping pushes the Fermi energy far 
from the Dirac point. For each position in space, one can 
map out the energy of the Dirac point E^. Shifts of E^, 
away from the average potential V gives the screened dis- 
order potential landscape whose autocorrelation function 
is C(r) = ((V(r) - V)(V(0) - V)). Theoretical calcula- 
tions for the correlation function are shown in Fig. |4] for 
different (average) electron carrier densities n using the 
Thomas-Fermi approximation and solving for the total 
density n to t — n + n e g self-consistently. The effect of 
electron doping is to slightly reduce both the magnitude 
of the potential fluctuations C(0) as well as the spatial 
correlation length of the puddles. 

Our predictions for screening and carrier transport in 
disordered 2D TI surface states should be testable in fu- 
ture experiments. In particular, the band asymmetry 
gives rise to interesting and qualitatively novel effects, 
e.g., causing the Dirac point, the charge neutrality point, 
and the minimum conductivity to occur at different car- 
rier densities. In addition, electron and hole 2D transport 
in TIs should manifest strong asymmetry, and in some 
situations with strong disorder (i.e., large rii mp or small 
d) there may not be any minimum conductivity plateau 
associated with the Dirac point. 

This work is supported by US-ONR and NRI-SWAN. 
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